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On the Hadronic Contribution to Light-by-light Scattering in — 2. 
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ABSTRACT 



We comment on the theoretical uncertainties involved in estimating the hadronic effects 
on the light-by-light scattering contribution to the anomalous magnetic moment of the 
muon, especially based on the analysis and results of Ref. 1. From the point of view of 
an effective field theory and chiral perturbation theory, we suggest that the charged pion 
contribution may be better determined than has been appreciated. However, the neutral 
pion contribution needs greater theoretical insight before its magnitude can be reliably 
estimated. 



The construction of an extremely high-precision experiment to determine the anoma- 
lous magnetic moment of the muon, a M = {g p — 2)/2, is underway at Brookhaven National 
Laboratory (BNL).* The anticipated design sensitivity, Aa p = 4 x 10~ 10 , will, if achieved, 
be about 20 times better than the results from CERN, corresponding to a magnitude 5 
times smaller than the Standard Model weak corrections. If the theory of known contri- 
butions is sufficiently good, such a measurement would either constrain or reveal physics 

[3 — 5] 

beyond the Standard Model. To draw such an inference, or even to check the weak 
radiative correction, requires a highly accurate determination of hadronic contributions. 
The 0(a 2 ) correction to a p coming from the O(a) contribution to hadronic vacuum po- 

[31 

larization to the photon propagator, has received a great deal of attention since it is 
approximately 35 times larger than the weak correction. The systematic error on this 
contribution can be reduced to the level of accuracy of the BNL experiment by a more 
precise low-energy measurement of R = a(hadron) jo{yT at VEPP-2M, w currently 
operating, or, in the future, at DA$NE. l6] 

There are other hadronic contributions in 0(ct 3 ) coming from the 0(a 2 ) contribution 
to the four-photon vacuum polarization tensor IL up \ a , the so-called contribution from light- 
by-light scattering. At first sight, since hadronic interactions are strong, this would appear 
to be very difficult to determine accurately, but Kinoshita et.al. [1] have estimated this 
contribution to be 49 x 10~ n with an quoted accuracy of about 10%. If correct, this 
contribution is a bit larger than the anticipated accuracy of the experiment, but four times 
smaller than the Standard Model weak correction. The method of calculation employed 
elementary pion contributions to the hadronic amplitude, together with vector meson 
dominance (VMD) of the photon couplings to the pion. They noted that their results 
happened to agree with a one-loop calculation involving elementary quarks with constituent 
quark masses, but they did not place great store by this second method because it is 
sensitive to the choice of the quark mass and would not appear to be a correct physical 
picture for momenta on the order of m M or so. 

Barbieri and Remiddi have recently raised doubts about the degree of certainty to be 
associated with the magnitude of the hadronic contribution in light-by-light scattering. As 
this is important for the interpretation of the experiment and crucial for drawing inferences 
about potential new physics, it is the purpose of this note to stimulate further discussion 

* For a summary and review of earlier measurements, see, e.g., Ref. 2. 
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of these important issues. 

The criticism in Ref. 7 is two- fold: (1) It is argued that the estimation in Ref. 1 by a 
quark loop is untrustworthy since it is so sensitive to the value chosen for the (constituent) 
quark mass. In this respect, they are actually in agreement with Ref. 1, who did not base 
their errors on the agreement with the estimation via a quark loop. So this should not be 
a point of dispute. (2) They doubt the accuracy claimed in Ref. 1, because it results "from 
the cancellation of different (gauge invariant) contributions, each of which is larger than 
the electroweak correction itself, In this regard, we believe that they are mistaken 

in suspecting that the approximation is unstable. In fact, the separate contributions, 
labelled A, B, and C in Ref. 1, while gauge-invariant with respect to the external photon, 
are not gauge-invariant contributions to U up \ a ; in particular, Bose symmetry among the 
four photons, while true of the sum, is not respected by the three separate sets of diagrams. 
Although each set, A, B, and C, yields a finite contribution to a p , it commonly happens in 
the calculation of radiative corrections that there are large cancellations among classes of 
diagrams that are not separately gauge invariant, and such an occurrence in and of itself 
is not necessarily cause to distrust the final result. 

Let us consider, ab initio how one may approximate the hadronic contributions to 
n^pAcr- The relevant momentum scale for the external photon momenta attached to H up \ a 
is set by the muon mass, an energy scale small compared to typical hadronic scales, with 
the exception of the pion. Therefore, it should be a good approximation to describe this 
by an effective field theory involving pions and photons only. At low energies, the pions 
may be regarded as the Goldstone bosons associated with chiral symmetry breaking, whose 
self-interactions are described by the 0(3) non-linear sigma model, gauged with respect to 
the Ui of electromagnetism. As a result, the interactions among pions at low-energy are 
actually weak, the scale of the chiral perturbation expansion being set by the mass m p of the 
p meson or thereabouts, well above the mass of the muon or pion. As a first approximation, 
then, one may neglect the self-interactions of pions entirely. For the charged pion, this is 
precisely what was done in Ref. 1, but one must understand that this is not a model for 
U up \ a , but the first term of a systematic expansion. The result for this contribution in 
Ref. 1 was given in Eq. (3.20). To illustrate that there should be no mystery here, a back- 
of-the-envelope estimate for this contribution to a p is as follows: all such contributions 
are proportional to m 2 : one power comes from the definition of a p ; the other comes from 
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the observation that the photon couplings are all chiral conserving whereas the anomalous 
moment requires a chiral-flip, so a mass insertion is required. This contribution will then 
be proportional to {m p /m^)^ 2 . There are 6 factors of the coupling e and a factor of 1/16/t 2 
for each loop. While each graph contributes a logarithmic subdivergence to Ii vp \ a , we know 
from gauge invariance that this divergence is cancelled in the sum. So we estimate 

(m p /m n ) 2 (a/47rf w .01(a/yr) 3 . (01) 

The actual result of — .04(a/"7r) 3 , given in Eq. (3.20) of Ref. 1, supports the view that, 
even though there are 21 graphs to be added,* the cancellation among divergent pieces 
also holds for the finite remainder. Thus, the final result is no bigger than our estimate, 
reinforcing our belief that the cancellation between subsets of graphs is simply a gauge 
cancellation. We do not have an argument for the sign of the contribution. 

Leaving aside the matter of the chiral anomaly for a moment, what are the corrections 
from interactions among the pions? The pion interactions do not become strong until 
one reaches to vicinity of the p-meson, so one would expect such contributions to g p — 2 
to be suppressed compared to the free pion contribution by a factor of approximately 
m 2 /m 2 4%. In Ref. 1, some of these corrections are introduced by VMD for the photon 
propagators, and it was found that they make corrections of more than 30% to various 
groups of diagrams contributing to a p and reduce the final answer for the charged pion 
contribution by a factor of more than 3. This casts a seemingly reliable approximation 
scheme in doubt. The VMD approximation used in Ref. 1 simply multiplies H up \ a by 
a common factor. Therefore, U up \ a remains gauge invariant,^ even th ough this VMD 
prescription does not respect the Ward identities for the couplings of photons to charged 
pions when off-mass-shell. Since H up \ a produces a convergent integral for the free pion 
loop, it is hard to understand how the VMD approximation, modifying the integrand on 
scales 1 arge compared to the pion mass, can change the final answer by a factor of three. 
This is a puzzle which suggests that the numerical integration be confirmed. It would add 
great confidence to the error estimates if the VMD corrections could be recalcul ated and 

* Actually there are only 8 different integrals; see Fig. 5 of Ref. 1. 
t I thank T. Kinoshita for emphasizing this point. 

* In fact, chiral perturbation theory has been extended to inclu de the p, but it is not clear how that 
formalism would be helpful here. 
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found to be of order rn%./m 2 p as expected. Of course, the integrand is not positive definite, 
so if the low-momentum region contributed little to the final answer, then modification of 
the high-momentum regime conceivably could produc e a large percentage change in the 
final result. However, in that case, it may be important to have a gauge-invariant VMD 
model. For the time being, we tend to trust the free pion result more. 

Another approach to estimating corrections would be to determine the effects of the 
higher order terms in the chiral perturbation expansion. That would be a complicated 
calculation that we have not attempted. However, we may anticipate one feature of these 
corrections, viz., because they involve higher-derivative interactions, they will lead to di- 
vergent contributions to a^. That simply signifies that, when the muon is included in the 
effective field theory, one must include a "direct" contribution to g p — 2 for which such a 
divergence provides a renormalization.* That is to say, the effective field theory contains 
a term of the form 

(02) 

where ay represents some effective coupling constant and A represents the scale at which 
this effective field theory breaks down. In the present context, we would expect A to be on 
the order of the masses of the vector mesons m p and m w . It has been tacitly assumed in 
Ref. 1 that the direct contribution is negligibly small. ^ Strictly speaking, one requires a 
renormalizable description of strong interactions in order to be able to assume that the bare 
ad is zero. QCD is, of course, such a theory but, unfortunately, we cannot reliably calculate 
its behavior in the low-energy regime relevant here. Even if we could, one must expect 
the physics omitted from any approximate description, renormalizable or not, to produce 
a direct interaction of this sort. For example, the weak correction itself may be regarded 
as a direct contribution at scales below M\y. [9] One might explore various renormalizable 
models for chiral-symmetry-breaking that are somewhat simpler than QCD; for example, 
one might choose the linear sigma model, gauged with respect to electromagnetism. This 
effectively provides a cutoff at the mass of the cr, and the sensitivity of the result to a 
mass (or the strength of the pion self-coupling) may provide an indication of the size of the 

* This sort of contribution to a M and the corresponding effective field theory is described in detail in 

Ref. 5. It occurs also for the tt° contribution; see below, 
t We shall revisit this assumption below in connection with the neutral pion contribution. 
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uncertainty in the corrections; however, extracting the contribution to a p from a two-loop 
contribution to IL up \ a is a bit daunting. 

In sum, we suggest that the result given for the non-interacting pion loop, Eq. (3.20) 
of Ref. 1, may be a reliable estimate of the charged pion contribution to an accuracy of 
about 4%. The VMD model is not gauge-invariant, although it does leave H up \ a gauge 
invariant. As a result, we would be inclined to trust the non-interacting result more than 
the somewhat smaller one favored in Ref . 1 resulting from their implementation of VMD. 
We now turn to consider the other contribution at the pion mass scale, the contribution 
of the 7T°. 

The effective field theory must be amended to include the effect of the chiral anomaly 
by adding 1101 

G-^-tPf^ (03) 

in the linear representation or the corresponding expression for the nonlinear realization 
in chiral perturbation theory.' 111 (Here, the pion decay constant f n ~ 93 MeV.) Eq. (03) 
describes the coupling of the 7r° to two photons at low energy, and this local vertex was 
used as the first approximation in Ref. 1 to the contribution of the 7r° to I\. vp \ a . The 
coupling constant G is predicted to be 1 in QCD in the "chiral limit" in which the pion 
is massless, and the experimental value agrees extremely well with this, to about 0.25%. 
Power counting shows that this local approximation to the interaction leads to a (loga- 
rithmically) divergent result for a p . This point-like interaction will be damped on scales 
associated with pion substructure, so, at scales on the order of m p , this local approxima- 
tion to the coupling G must be modified. This is qualitatively equivalent to introducing 
damping factors via VMD as in Ref. 1, with a result that is given in their Eq. (3.27). In 
this case, their approximation to the vertex is gauge-invariant, and, therefore, this can be 
expected to provide a believable estimate of the size of this correction. The result should 
not be very different from simply introducing a cutoff on the divergence at the scale m p , 
and, since the divergence is only logarithmic, the result will not be very sensitive to the 
precise value chosen. Again, a back-of-the-envelope calculation gives for any one diagram 

~ <^ 2e2 ^ 2 Hrn 2 p /ml) = \(£f) 2 ^) 3 Hrf/ml) - .006(a/ 7 r) 3 . 
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There are two equal pairs of two diagrams each,* so we ought to multiply this by at least 
a factor of two; that is still about a factor of 4 less than the result of .05(a//r) 3 given in 
Eq. (3.27) of Ref. 1, using VMD to cut off the divergence. 

It is actually quite surprising that the magnitude of the 7r° contribution is as large as 
the charged pion contribution.* Even though it is obviously the same order in a, it is in 
effect one- loop order higher. That is, the anomalous coupling a/nf^ itself is a one-loop 
contribution, and it enters squared in a two-loop graph. Thus, this is in effect a four- 
loop contribution, to be compared with the three-loop contributions involving the charged 
pion. So, a priori, one would anticipate this contribution being suppressed by 1/167T 2 . Even 
though there are many more diagrams in the charged pion case, the cancellations among 
diagrams as required by gauge invariance arranges for the final result to be no bigger than 
our estimate. Thus, whereas the charged pion contribution is somewhat smaller than one 
might have guessed given the number of diagrams, the 7r° contribution is somewhat larger. 
This enhancement may be attributed to a host of small factors: the anomaly is twice as 
large as we would have guessed and it enters squared; the log enhancement is about a factor 
of 3; a factor of 2 from (m n / f n ) 2 ; finally the various diagrams may add constructively in 
this case and destructively in the charged pion case. Thus, remarkably, the additional 
factor of 1/167T 2 may be largely compensated, but it remains rather surprising and would 
be worth double-checking. 

Combining the neutral pion result with the charged pion loop contribution (Eq. (3.20) 
rather than Eq. (3.24) of Ref. 1,) one obtains the result 

(04) 

coincidentally about 3 times smaller than that given there. This depends upon the signif- 
icant cancellation between the charged and neutral pion contributions. 

If this were the whole story, the situation would be quite satisfactory. However, the 
fact that the local coupling Eq. (03) results in a ultraviolet divergence implies that there 
must also be a "direct" contribution to a M coming from physics above the cutoff (m p ) 
whose natural size cannot be much smaller than the radiative correction due to the 7r°. 

* The two diagrams of Fig. 5 of Ref. 1 must be added to two having the muon line reversed. 

* According to Ref. 1, the free pion result is about equal but of opposite sign. 



a M (had2) = 0.014 



a 



IT. 
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There are several ways to see this: If we regard as a bare coupling, then the divergence 
of the 7T° contribution renormalizes it. Alternatively but equivalently, we may infer from 
the logarithmic divergence a contribution to the (3 function for the renormalized coupling 
ad(M), where M is the renormalization scale. It would be unnatural, in the technical sense 

[12] 

of the word, for it to be small compared to the size implied by varying M a bit over the 
range of interest. Since h^mp/m^) ~ 1.6, we would expect the direct term to be not much 
smaller than the 7r° contribution calculated in Ref. 1. This suggests that important physics 
has been omitted from this method of calculation whose magnitude may well be comparable 
to the contribution calculated. 

Of course, we expect that the 7r°77 vertex, even for the simple triangle diagram, will 
be damped whenever one of the external momenta becomes large compared to the typical 
scale of strong interactions. On the other hand, the Adler-Bardeen theorem' 101 guarantees 
that, for zero external momenta, the strong interactions do not in fact modify the anomaly 
associated with the underlying fundamental quark-loop diagram. The theoretical challenge 
is to figure out how to move away from the local approximation in a manner that is 
sufficiently well- controlled to provide an estimate of the error on the calculated contribution 
to g^ - 2. 

The implications of this uncertainty associated with the 7r° contribution are severe: If 
the data turn out not to agree with the theoretically predicted value and the experimental 
errors are sufficiently small, then, rather than necessarily a signature of new physics, the 
deviations could just as well be ascribed to an erroneous assumption about the magnitude 
of the hadronic contribution from light-by-light scattering. In this respect, we agree with 
Ref. 7 that more work needs to be done. 

ACKNOWLEDGEMENTS 

I have enjoyed discussions on this subject with T. Kinoshita and D.R.T. Jones. This 
work has been supported in part by the U.S. Department of Energy. 



7 



REFERENCES 



1. T. Kinoshita, B. Nizic, and Y. Okamoto, Phys. Rev. D31, 2108 (1985). 

2. F.J.M. Farley and E. Picasso, "The Muon g — 2 Experiments" in Quantum Electro- 
dynamics ed. T. Kinoshita, Singapore: World Scientific, 1990. 

3. T. Kinoshita and W.J. Marciano, "Theory of the Muon Anomalous Magnetic Mo- 
ment" in Quantum Electrodynamics, op.cit. 

4. T. Kinoshita, Talk at 10th International Symposium on "High Energy Spin Physics," 
Nagoya, Nov. 9-14, 1992, (CLNS 93/1186.) 

5. C. Arzt, M.B. Einhorn, and J. Wudka, Michigan UM-TH-92-17 (Phys. Rev. D, to be 
published.) 

6. The Da$ne Physics Handbook ed. L. Maiani et.al., Frascati, INFN, 1992. 

7. R. Barbieri and E. Remiddi, The DA&NE Physics Handbook, op.cit. vol. II, p. 301. 

8. G. Ecker et.al., Phys. Lett. B223, 425 (1989). 

9. E. Witten, Nucl. Phys. B122, 109 (1977). 

10. S.L. Adler and W.A. Bardeen, Phys. Rev. 182 1517 (1969). 

11. J. Wess and B. Zumino, Phys. Lett. 37B, 95 (1971), E. Witten, Nucl. Phys. B223 
422 (1983). 

12. S. Weinberg, Physica 96A 327 (1979). 



8 



